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Abstract 

Inclusive branching ratios involving b — > r transitions are measured in ap- 
proximately four million hadronic Z decays collected by the ALEPH detector at 
LEP. The fully- inclusive branching ratio BR(b — > T~i>T-X) and the semi-inclusive 
branching ratio BR(b are measured to be (2.43ib0.20ib0.25)% and 

(0.88ib0.31ib0.28)%, in agreement with the standard model predictions. Upper 
limits on the branching fractions BR(B^ t^i^t) and BR(b svu) are set to 
8.3 X 10"'' and 6.4 x lO-"^ at the 90% C.L. These results allow a 90% C.L. lower 
limit of 0.40 (GeV/c^)~^ to be set on the tan/3/mjj± ratio, in the framework of 
type-II two-Higgs-doublet models. 
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1 Introduction 



Third-generation fermions are involved in both the initial and final states of b — > r 
transitions. A measurement of branching ratios pertaining to these transitions can 
be compared to the standard model predictions, yielding direct tests of heavy fermion 
interactions. These interactions are especially sensitive to the mechanism underlying 
the electroweak symmetry breaking, i.e., to the origin of mass, and allow extensions of 
the standard model to be constrained in this respect. In the study presented here, the 
branching fractions of the following processes are considered. 

• The inclusive branching fraction BR(b r^z/^X) can be compared to the standard 
model prediction of 2.30ib0.25% as computed in the framework of the Heavy 
Quark Effective Theory (HQET). The standard model transition, illustrated in 
Fig. [l|a, could be modified by the exchange of a new charged boson, as shown in 
Fig. |I]b with a charged Higgs boson. In any type-II two-Higgs-doublet model, such 
as the minimal supersymmetric extension of the standard model, the corresponding 
contribution to the transition amplitude is proportional to {tan (3 / m^±)'^ ||, ^, ||, 
where tan /3 is the ratio of the vacuum expectation values of the two Higgs doublets, 
and mpji is the mass of the charged Higgs boson. 

• A similar test can be performed with the semi-inclusive decay b t~ i?T-D*~^X. In 
the standard model, this branching fraction is expected to be approximately 1% 
The proportion of D* relative to D final states would be reduced in the presence of 
a charged Higgs boson [||. 

• The exclusive decay mode B~ — t~Uj. (Fig. |l|c) has a branching fraction predicted 
to be 7.4 X 10-5 (/B/160MeV)^ (|Vub|/0.004f in the standard model §. In type-II 
two-Higgs-doublet models (Fig. |]d), the rate is enhanced by the multiplicative 
factor [m|_ (tan/3/mH±)^ — 1]^. 

• Although it does not involve a b — r transition, the decay b — > suu, with a bran- 
ching fraction predicted to be of the order of 5x IQ-^ in the standard model (Fig. |l|e), 
can also be exploited to put constraints on a variety of extensions of the standard 
model (Fig. |l]f) as advocated in Refs. 

The measurements of the branching fractions of these four processes, b r^z/^X, 
b r^pT-D*^X, —>■ T^Vj. and b suu, are performed with all data collected by the 
ALEPH detector between 1991 and 1995, at centre-of-mass energies at and around the Z 
resonance, corresponding to approximately four million hadronic Z decays. As described 
in previous ALEPH studies applied to smaller data samples |10|, |ll|, a single technique is 
used to extract the four branching ratios. The initial-state b quark is identified by means 
of standard b-tagging algorithms, and the final-state r lepton (or the two neutrinos in the 
case of b — > suu) is identified by the missing energy carried away by the neutrinos. 

The energy-fiow and b-tagging algorithms are briefiy described in Section ^ along 
with properties of the ALEPH detector relevant for the present study. The selection 
algorithms aimed at rejecting background events with a large missing energy are presented 
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in Section |^. The method to estimate the residual background yields is detailed in 
Section^. The branching ratios BR(b t^z/t-X) and BR(b t^Ut-D*^X) are determined 
in Section ^, and the limits on BR(B~ —>■ T~j?r) and BR(b — > suu) are extracted in 
Section ^ Finally, an alternative measurement of BR(b —>■ T~i?r^) with opposite-sign 
di-leptons in the final state is presented in Section 0. The results are interpreted in the 
framework of type-II two-Higgs-doublet models in Section ^ and summarized in Section ^ 




a) 



H- / 



b) 
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Figure 1: Feynman graphs for the processes b ^ r i^rX, B t Vr and b — > svv^ in the standard 
model (a,c,e) and in two-Higgs-doublet models with the exchange of a charged Higgs boson (b,d,f). 



2 The ALEPH Detector 



A detailed description of the ALEPH detector can be found in Ref. [0], and of its 
performance in Ref. [|13|. Charged particles are detected in the central part, consisting 
of a precision silicon vertex detector, a cylindrical drift chamber and a large time 
projection chamber, measuring altogether up to 31 space points along the charged particle 
trajectories. A 1.5 T axial magnetic field is provided by a superconducting solenoidal 
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coil. Charged particle transverse momenta are reconstructed with a. l/px resolution of 
(6 X 10~^ 5 X \{)~^ /px) (GeV/ c)~^ . In the following, good tracks are defined as charged 
particle tracks reconstructed with at least four hits in the time projection chamber, 
originating from within a cylinder of length 20 cm and radius 2 cm coaxial with the beam 
and centred at the nominal collision point, and with a polar angle with respect to the 
beam such that I cos 01 < 0.95. 



Jets originating from b quarks are identified with a lifetime b-tagging algorithm [|T4|, 
which takes advantage of the three-dimensional impact parameter resolution of charged 
particle tracks. For tracks with two space points in the silicon vertex detector {i.e., 

1 cos 6' I < 0.7), this resolution can be parametrized as (25 + 95/p) /im, with p in GeV/c. 

In addition to its role as a tracking device, the time projection chamber also measures 
the specific energy loss by ionization dE/dx. It allows low momentum electrons to be 
separated from other charged particle species by more than three standard deviations. 

Electrons (and photons) are also identified by the characteristic longitudinal and 
transverse developments of the associated showers in the electromagnetic calorimeter, 
a 22 radiation length thick sandwich of lead planes and proportional wire chambers with 
fine read-out segmentation. The relative energy resolution achieved is 0.18/V^ {E in 
GeV) for isolated electrons and photons. 

Photon conversions to e"'"e~ in the detector material are identified as a pair of 
oppositely charged identified electrons satisfying the following conditions: (i) the distance 
of closest approach to the beam of the two reconstructed tracks is larger than 2 mm; (ii) 
the distance between the two tracks at their point of closest approach is smaller than 

2 cm in space; (in) the position of this point is consistent with a material boundary; and 
(iv) the invariant mass is smaller than 20 MeV/c^, when calculated as for an e"'"e~ pair 
coming from this point of closest approach. 

Muons are identified by their characteristic penetration pattern in the hadron 
calorimeter, a 1.5 m thick yoke interleaved with 23 layers of streamer tubes, together with 
two surrounding double-layers of muon chambers. In association with the electromagnetic 
calorimeter, the hadron calorimeter also provides a measurement of the hadronic energy 
with a relative resolution of 0.85/-\/£' {E in GeV). 

Taus are identified by the missing energy carried away by their decay neutrinos. The 
total visible energy is measured with an energy-fiow reconstruction algorithm which 
combines all the above measurements, supplemented by the energy detected at low 
polar angle (down to 24 mrad from the beam axis) by two additional electromagnetic 
calorimeters, used for the luminosity determination. The relative resolution on the total 
visible energy varies between 0.60/v^ for high multiplicity final states and 0.25/v^ 
for final states of low multiplicity without neutral hadrons. In addition to the total 
visible-energy measurement, the energy- flow reconstruction algorithm also provides a list 
of reconstructed objects, classified as charged particles, photons and neutral hadrons, and 
called energy-flow particles in the following. 
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3 Event selections 



A b ^ r transition followed by the decay of the r always produces two energetic i^t-'s. 
Such transitions can therefore be identified in e"*"e~ bb events on the basis of a large 
measured missing energy. The argument holds as well for b —>■ suu decays, which can be 
selected along the same lines as b ^ r transitions. In this section, the algorithms aimed 
at rejecting background events leading to large measured missing energy are described. 
An alternative selection, based on the presence of two identified leptons with opposite 
electric charge, is described in Section ^ 



3.1 Preselection 

Approximately four million hadronic Z decays are selected in the data collected between 
1991 and 1995 at energies at and around the Z resonance when the standard criteria (at 
least five good tracks carrying at least 10% of the centre-of-mass energy [1^) are applied. 
To keep only two-jet events well contained in the detector acceptance, the polar angle is 
required to satisfy | cos6'thrust| < 0.7 to match the acceptance of the vertex detector, and 
the thrust value must exceed 0.85. 

Each event is then divided in two hemispheres with respect to the plane perpendi- 

1 2 

cular to the thrust axis. In each hemisphere, the missing energy -E^iss is defined as the 

12 12 

difference between the expected true energy -EtrW and the measured visible energy E^l^. 
The latter is determined from the total energy of all energy-flow particles contained in 
that hemisphere, while the former is estimated from the centre-of-mass energy ^/s and 
with energy-momentum conservation: -E^'ue = (s -|- 2 ~ "^2,i)/2a/s, where mi and m2 
are the measured invariant masses in the two hemispheres p!6| . 

The main background to the final states searched for in this analysis consists of Z — > bb 
or cc events followed by a semi-leptonic b or c decay into an electron or a muon (hereafter 
called lepton and denoted £), with a large missing energy carried away by the neutrino 
ui. This background can be considerably reduced by rejecting hemispheres in which 



a lepton is identified. The standard lepton identification |T3] is not used here, but is 
replaced by much looser criteria. A good track is tagged as an electron either if the 
dE/dx is compatible with that of an electron and incompatible with that of a pion, or 
if the transverse and longitudinal profiles of the associated electromagnetic shower are 
compatible with those of an electron. Electrons and positrons originating from identified 
photon conversions are not considered in this process. Similarly, a good track is tagged as 
a muon if it is associated to a few hits either in the last layers of the hadron calorimeter 
or in the two layers of muon chambers. A lepton identification efficiency over 95% is 
achieved for electrons (muons) with a momentum in excess of 1 (2) GeV/c. 

In addition, to reduce the contamination from Z decays into lighter quark pairs (uu, dd, 
ss and cc) in which a large missing energy is faked due to finite detector resolution effects, 
the final sample is enriched in Z ^ bb events with b tagging. The corresponding criteria 
are dependent on the signal final state and are described in the following subsections. 
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3.2 The b T-i?^^ final state 

For the b t~iJt-X final state selection, a hemisphere is kept only if the opposite 
hemisphere is tagged as arising from a b quark. Each good track in the opposite 
hemisphere is assigned a probability of originating from the primary interaction point, 
on the basis of its impact parameter significance. The confidence level a'^'^™ that all good 
tracks in that hemisphere come from the primary interaction point, determined under 



the assumption that the individual probabilities are uncorrelated |T^, is required to be 
smaller than 1%. 

Residual backgrounds like Z — > r+r^ decays, two-photon processes or beam-gas 
interactions are likely to yield a large missing energy in the final state and might therefore 
bias the analysis. These background sources are reduced down to negligible levels by 
requiring at least seven good tracks and a total missing energy smaller than 50 GeV, with 
almost no additional loss of the signal. 



3.3 The B — ¥ r Ut and b — > suv final states 

The same b-tagging technique as for the b r^z/^-X selection is applied but the residual 
backgrounds are rejected in a slightly different manner. The opposite hemisphere is 
required to contain at least six good tracks, and its missing energy must be smaller than 
25 GeV. This latter cut reduces the effect on the missing-energy distribution in the signal 
hemisphere, which is expected to extend towards larger values than in the previous case. 

3.4 The b r-i/^D*±X final state 

The b r^z/T-D*^X final state selection is primarily based on an exclusive reconstruction 
of a D*^, very similar to that performed in Ref. |TB[ although with slightly tighter cuts 
as a result of the larger statistics accumulated until 1995. 

Candidate D*^'s are searched for in the decay channel D*^ — D*^7rsoft, followed by 
D° ^ K-7r+7r-7r+, D° ^ K-7r+, D° ^ K-tt+tt" or D° ^ Kgyr+vr-, selected as follows. 

Neutral pions are identified as photon pairs with an invariant mass compatible with 
the 7r° mass and a probability in excess of 1%. Candidate Kg's are identified as pairs 
of oppositely charged particle tracks (assumed to be charged pions) forming a secondary 
vertex at least 2 cm away from the primary interaction point, and with a reconstructed 
invariant mass within 5MeV/c2 of the Kg mass. Candidate K^'s are identified as good 
tracks of momentum in excess of 3GeV/c and with a dE/dx compatible with that of a 
charged kaon within two standard deviations. In addition the angle 9* between the kaon 
momentum direction, evaluated in the D*^ rest frame, and the D° line-of-flight must satisfy 
I cos r I < 0.9. 

The reconstructed invariant mass of the candidate D°'s is required to be within 
15MeV/c^ (50MeV/c^) of the D*^ mass for the first (last) two decay channels. The ratio 
of the D° energy to the beam energy is required to be between 0.25 and 0.50. In addition. 
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the reconstructed D° vertex must be incompatible with the primary interaction point by 
at least two standard deviations. 

Each reconstructed D° is associated with a low momentum pion {Ptt^^i^ < 4.2GeV/c) 
to form a D*^. The reconstructed invariant mass difference m£)*± — moo is required to be 
within 20MeV/c^ of the nominal mass difference. 

Finally, the D*^ hemispheres are tagged as arising from a b quark by requiring the 
presence of an additional good track in a cone of half-angle 30° around the D° momentum 
direction. This track must be incompatible with originating from the primary interaction 
point by more than three standard deviations, and its electric charge must be opposite to 
that of the vTsoft- No b tagging needs therefore be applied to the opposite hemisphere. 



3.5 Event yields 



The expected event yields were determined with the following samples of events, simulated 



with JETSET [^, and subsequently processed through the full GEANT simulation 
of the detector: 



approximately four million hadronic Z decays for the background estimate; 

over 40 000 bb events, each containing a b ^ t^z/t-X, with the r polarization 
determined in the limit of the free-quark model |l|] for mb = 4.8GeV/c^ and 
mc = 1.4GeV/c2; 

about 20 000 bb events with at least one r^v^- and 100% polarized r's ; 

almost 20 000 bb events with at least one b — > sz/i/, generated with the B — > XgZ/i/ 
decay probability given in Ref. [^, in which the new-physics parameters were 
conservatively chosen to minimize the expected selection efficiency. The model was 
improved with a realistic Xg invariant mass distribution, described as a Gaussian of 
mean and variance 1.35 and 0.60GeV/c^ 0, supplemented by two peaks at the K 



and K* masses 22| with branching fractions of 0.06 ± 0.02 and 0.29 ± 0.01. 



In all these samples, the polarization of the r's was properly taken into account in the 
simulation of their decay kinematics and events containing a Ab were reweighted 
according to the measured Ab polarization [^, so that the correct distribution of the Vj- 
momenta, and therefore the missing energy, be obtained. 

The numbers of hemispheres expected to be selected from the various background and 
signal processes and the numbers of hemispheres selected in the data are listed in Table |l[ 

The measurement of the various branching fractions relies on the observation of an 
excess of events at large -Emiss in the missing-energy distribution of the data with respect 
to that of the simulated background. Given the limited number of events expected from 
the various signal processes (Table |l|), it is essential to have a detailed understanding of 
all uncertainties affecting the background distribution both in shape and normalization 
(Section ||) and of the signal selection efficiencies (Section ||). 
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Table 1: Numbers of hemispheres selected in the data in the four different final states. Also indicated 
are the numbers of hemispheres expected from the various background and signal processes. For the 
latter, the standard model branching ratio values were assumed. The standard model branching ratios 
and the selection efficiencies are indicated in brackets. 



Final state selection 




r-z/^D*±X 


r u^-, suu 


Data 


Hemispheres selected 


166342 


1464 


156910 


Simulation 


Hemispheres expected 


162456 


1538 


153093 


Background (details) 


u, d, s 


9215 


1172 


8855 


b, c without leptonic decay 


138118 


80 


129886 


b, c with leptonic decay 


10645 


180 


10091 


Signal (details) 


b ^ r-z/,X (2.3%) 


(12%) 4478 




4247 


b ^ r-z/,D*±X (1.0%) 




(0.34%) 81 




b t~Ut-X without D*^ 




25 




B- ^ T'Vr (7.4 X 10^^) 






(8.1%) 4 


b suu (5.0 X 10"^) 






(8.8%) 10 



4 Background estimate 

The missing-energy distribution of the background to b — r and b suu transitions 
depends on detector performance in three areas, (i) the visible-energy reconstruction for 
hemispheres with a priori no missing energy such as Z hadronic decays into light-quark 
pairs; (ii) the effectiveness of the lepton veto to reject background from b, c — >• iuiX; 
and (Hi) the performance of the b-tagging. To minimize the influence of possible 
inaccuracies in the detector simulation, these three quantities are derived from the data 
and subsequently used to re-weight the simulated events. 



4.1 The visible-energy calibration 

As mentioned above, the accuracy of the visible-energy reconstruction can be studied 
with a sample of events in which no significant missing energy is expected. Such a sample 



is selected in a way identical to that followed in Sections |3.1| and |3.2| except that the b- 
tagging criterion is replaced by a b-rejection one requiring a'^''™ > 0.5. It yields fractions 
of 88.7%, 10.5% and 0.8% for (uu,dd,ss), cc and bb events, respectively. 

The visible energy E^is is the sum of the energy Echa. of the charged particles 
reconstructed in the central tracker, the energy Ep^o of the photons detected in 
the electromagnetic calorimeter, and the energy E'neu of the neutral hadrons, the 
reconstruction of which makes use of all identification and energy measurement capa- 
bilities of the detector (Section 0). These three contributions to the visible energy 
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may therefore be affected by different systematic uncertainties and must be calibrated 
separately. The hemispheres were thus binned in a grid according to the fractional 
contributions to the visible energy of the charged, photon and neutral hadron components 
{Echa/E^is, -Epho/-Evis, E^en/E^is). The accuracy of the detector simulation was examined 
by comparing the observed and simulated missing energy in each of the bins. 

A noticeable inaccuracy of the simulation was found in hemispheres with a large 
proportion of neutral hadronic energy. This disagreement with the data could lead 
to a large systematic uncertainty in the final result. Moreover, the neutral hadronic 
energy is expected to be smaller in bb events than in other hadronic Z decays. It was 
therefore decided to reject all hemispheres for which E^eu is in excess of 7 GeV, both in this 
calibration procedure and in the selections described in the previous section, preserving 
69.6% of the events in the data and 68.7% in the simulation. 

The residual differences between the data and the simulation were corrected by scaling 
Echa., -Epho and -Enou by /cha, /pho and /neu, respectively, in the simulated events. The values 
of these calibration factors were obtained by minimizing the following x^, 

= J2{(^vis,i) - fcha{E^,a,i) + /pho(-E^o,i) + /neu (-£^^u,i) } / Crf , (1) 
i 

where the mean energy values are computed in each bin i, E^^^f is the total visible energy 
measured in that bin, the -E^x « values are the charged, photonic and neutral hadronic 
energies expected in that bin, and ai is the uncertainty due to the limited statistics of 
the event samples. The result of the fit is shown in Table 0. While /cha and /pho are 
compatible with unity, a sizeable calibration is found to be still necessary for the neutral 



hadronic energy, despite the cut at 7 GeV. This effect, already reported in Ref. ||I0|, results 
from an inadequacy of the simulation of nuclear interactions. 



Table 2: Calibration parameters (see text for a detailed description) for Evh 



/cha 


1.002 ± 0.001 


/pho 


1.000 ± 0.002 


/nou 


0.936 ± 0.010 


/b 


0.999 ± 0.001 



As another consequence, the -Eneu resolution is better in the simulation than in the 
data. In addition to the rescaling, the simulated neutral hadronic energy was therefore 
smeared on an event-by-event basis, leading to a (1.1 ± 1.6)% relative deterioration of 
the -Eneu resolution. This smearing was performed in order to equalize the observed and 
expected numbers of events with a reconstructed -Emiss in excess of 14 GeV. Indeed, in 
events where no real missing energy is expected, such a large reconstructed -Emiss value is 
often due to a loss of neutral hadronic energy. 

The distributions of the reconstructed missing energy in the control sample are 
displayed in Fig. H, both for data and simulated events, before and after the calibration 
procedure. The agreement is found to be adequate over the whole spectrum, and in 
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Figure 2: Reconstructed missing-energy distributions using the light-quark tag and e//i veto (a) before 
calibration; (b) after missing-energy recalibration and applying the cut i^nou < 7 GeV. 



particular at large missing energy where the background to the signal of interest in this 
paper has to be evaluated. 

Finally, due to the specific hadronization and fragmentation of b-quark jets, the sharing 
between the charged, photonic and neutral hadronic energy may lead to an overall visible- 
energy shift different from that observed for lighter quark species. The visible energy of 
hemispheres tagged as arising from a b quark was therefore rescaled by an overall factor 
/b, determined from the b-tagged sample of Section so as to match the peak positions 
(which are not affected by events with truly missing energy) of the expected and observed 
visible-energy distributions. This b-specific calibration factor is found to be compatible 
with unity (Table |^). 

4.2 The lepton rejection effectiveness 

In order to estimate the effectiveness of the lepton veto, aimed at rejecting Z — » bb or cc 
events followed by a semi-leptonic b or c decay into an electron or a muon, the veto was 
applied to unbiased, pure, lepton data and simulated samples. 

A 97.6% pure electron sample was obtained by selecting hadronic events with a photon 
conversion, as described in Section |^, in which only one of the two particles has to be 
identified as an electron. Photon conversions are well suited to study the electron veto 
since the hadronic environment is very similar to that of semi-leptonic heavy-fiavour 
decays. The electron-veto effectiveness is thus given by the probability, after background 
subtraction, that the second particle be identified as an electron by the lepton- veto criteria. 
However, this probability depends on (i) the electron momentum, the spectrum of which 
is different in photon conversions and in semi-leptonic b/c decays; and (ii) the number of 
wires used for the dE/dx measurement, which tends to be smaller in photon conversions 
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than in semi-leptonic b/c decays due to shared hits with the first electron of the converted 
pair. Therefore, the electron identification probability was mapped as a function of the 
electron momentum and the number of hits associated to the corresponding track, both 
for the data and the simulated events with photon conversions. The ratio of these two 
maps was then used to re- weight the simulated b/c events with semi-leptonic decays. 

For the muons, a sample originating from the 77 — ^ process was selected. The 

absence of hadronic environment is less important because the identification is achieved 
with the external layers of the hadron calorimeter and the muon chambers. A 97.1% 
pure sample of low-momentum muons was obtained by requiring exactly two good tracks 
(accompanied by no neutral energy-flow particle) with opposite electric charge, one of 
which being identified as a muon, a total momentum less than 30 GeV/c and an invariant 
mass smaller than 2GeV/c^. Above 3GeV/c, muons are expected to traverse entirely 
the hadron calorimeter and the muon chambers, making the identification efficiency 
independent of momentum. The 77 — fi'^^~ sample was therefore supplemented by 
Z decays into muon pairs to determine the latter. As above, and after background 
subtraction, the muon-veto effectiveness was given by the probability that the second 
particle be identified as a muon by the lepton-veto criteria as a function of momentum, 
both in the data and in the simulation. The ratio was then used to re-weight the simulated 
b/c events with semi-leptonic decays. 

The identification probability for electrons and muons, determined from these samples, 
is displayed as a function of the lepton momentum in Figs. ^ and §3, for data and 
simulation, and is compared to the identification efficiency for simulated b/c events with 
semi-leptonic decays. The data-to-simulation ratio, applied as a correction factor to the 
latter, is consistently found to be close to unity. 
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Figure 3: Lepton identification efficiencies (a) for electrons and (b) for muons as a function of momentum 
in the control samples (see text) for the data (open circles) and the simulation (full circles), and their ratio 
(triangles). The histograms show the identification efficiency for simulated b/c events with semi-leptonic 
decays. 
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4.3 The b-tagging efficiency 



The efficiency of the b-tagging criterion can be determined directly from the data using 



the double tag method [|1J], so as to minimize the systematic uncertainties related to the 
limited knowledge of b-hadron production and decay. The small correlation between the 
b-tagging probabilities in the two hemispheres is negligible for the present purpose. 

From the knowledge of the fractions i?b and Rc of hadronic Z decays into bb and 
cc, the efficiencies and of the b-tagging criterion on b and c hemispheres can be 
measured in the data by comparing the fraction /i of hemispheres which pass a given 
Q,hemi 1^^^ with the fraction /2 of events in which both hemispheres pass the same cut. 
These fractions are given by 

/i = -Rb^b + Rc^c + (1 — -Rc — -Rb)^x, 

/2 = R,,el + R,el + {l-R,-R^)el (2) 

where is the hemisphere tagging efficiency for light-quark hemispheres, the value of 
which is determined directly with the data as explained in Ref. [|I^. The difference 



between the observed and expected £x values for a given a'^^™' cut arises from from 
different impact parameter resolutions in the data and in the simulation. This difference 
was estimated from the fraction of hemispheres satisfying a^g™ < 0.01, where ctng™ 
is the hemisphere probability recomputed with only those tracks having negative impact 
parameter (and reversing the impact parameter sign). Such tracks nearly always originate 
from the primary vertex, and account for resolution differences. The relatively small 
contribution of long-lived particles (such as Kg or A°) is taken from the simulation. 

With the cut a'^'^™ < 0.01, was found to be 2.46%, i.e., slightly larger in the data 
than predicted by the simulation by a factor of 1.17 ± 0.11. With the value obtained 
this way, Eq. (^ was solved for e\, and Ec- The tagging efficiencies for b (56.7%) and c 
(14.8%) hemispheres were also found to be larger in the data by factors of 1.03 ±0.01 and 
1.03 ± 0.03, respectively, and were used to re-weight the simulated events accordingly. 



5 Measurement of b — > t z^t-X and b — > t i>'t-D*^X 



The b — > r~z/T-X and b — r~z/T-D* X branching fractions can be determined by evaluating 
the excess of events over the background expected at large missing energy from the 
two selections (Sections ^]2| and p.4|) . The -Emiss distributions in the b — > r^z/^X and 
b — s> r~z/T-D*^X final states are displayed in Figs. ^ and for the data and for the 
re- weighted simulation (Section ^), normalized to the number of observed events. 



5.1 Branching ratio extraction 



An excess of events at large -Emiss is indeed observed over the background in both 
distributions. The corresponding branching ratios were extracted with a binned-likelihood 
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Figure 4: Missing-energy distributions in the b t~Dt^ (a) and the b — > r~pT-D*='=X (b) final states, 
for the data (dots) and the simulation (histograms). The latter is subdivided in (i) the fitted signal 
contribution (empty histogram); (ii) the contribution from b and c semi-leptonic decays (light hatching); 
and (in) the residual background (dark hatching). Also indicated are the -Emiss intervals considered for 
the branching ratio measurements (Section 5.1). 



fit of tlie expected to tlie observed missing-energy distributions, keeping tfie normalization 
to tfie number of events observed. Tlie fit was performed in -Emiss intervals chosen so as 
to (i) minimize their total uncertainty, defined as the quadratic sum of the statistical and 
the systematic contributions (Section |5.2|) ; and (ii) make the measurements statistically 
independent of that of the r~z/^ branching fraction. The optimal intervals, 

indicated in Fig. ^ are found to be 18-35 GeV and 12-35 GeV, respectively. 

The numbers of events in these intervals, observed in the data and expected from 
both the background and the signal processes, are displayed in Table In particular, the 
cascade decay b —>■ D~X with D^T T^i?r, which yields an i^miss spectrum similar to that 
of the decay b r^z/^-X, is included in the backgrounds to this signature. 

These numbers yield measured branching fractions of 



BR(b^r-z/^X) = [2.43±0.20(stat.) ±0.25(syst.)]%, 

and 

BR(b r-z/^D*±X) = [0.88±0.31(stat.) ±0.28(syst.)]%. 



in agreement with the standard model predictions of (2.30 ± 0.25)% and approximately 
1%, respectively. The systematic uncertainties, also indicated above, are discussed in the 
next section. 
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Table 3: Numbers of events observed in the data at large E'miss in the b T~DrX and b t~ Dr^*^X 
final states. Also indicated are the events expected from the various background and signal processes. 
For the latter, the fitted branching ratio value was assumed. The corresponding selection efficiencies are 
indicated in brackets. 



Final state selection 




r-z/^D*±X 


Emiss interval 


[18,35] GeV 


[12,35] GeV 


Data 


Hemispheres selected 


2094 


162 


Expected background (details) 


u, d, s 


17 


2 


b, c with leptonic decay 


1001 


51 


b ^ D^^X, D3- ^ T~iJr 


214 




b, c without i, 


84 


36 


Expected signal (details) 


b ^ r Ut-X 


(2%) 778 




b ^ T-UrB*^X 




(0.26%) 63 


b T^Ur^ withoutD*=^ 
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5.2 Systematic uncertainties 

The measured branching ratios BR(b — r^z/^-X) and BR(b — ^ Ut-D*^X.) can be affected 
by several sources of systematic effects on the signal selection efficiencies and by residual 
inaccuracies of the background simulation. The relevant sources common to b ^ t^z/t-X 
and b —>■ Ur^*^X are addressed in turn below, and their effects on the branching 
fractions are summarized in Table ^ 



1. The missing-energy distributions of b — > r~z/^X and b, c — * iuiX. depend on the expe- 
rimentally determined quantities (xb), (xc), BR(b ii'iX) and BR(b c — > iuiX). 
However, the small, residual differences between the i^miss distributions in the 
data and the simulation, obtained by reversing the lepton-veto criteria (Fig. |^), 
were found to be well accounted for by the measured uncertainties on their 
values. The whole extraction procedure of the branching fractions was therefore 
repeated by varying these quantities within their uncertainties, so as to evaluate the 
corresponding systematic effects. Similarly, the branching ratio dependence on the 
D** content in b decays was accounted for according to Ref. . 



2. The b-fragmentation modelling affects the b-energy spectrum, and therefore also the 
missing-energy spectrum in b — >■ t~Ut-X and b, c — >■ ii'iX. The effect of its knowledge 
was estimated in a similar manner by using the four different fragmentation functions 
described in Ref. [|^, adjusted to reproduce the measured value of (xb). The 
corresponding systematic uncertainty was defined as the largest change in the 
measured branching fractions. 
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Figure 5: -Emiss distributions obtained for the data (dots) and the simulation (histogram) from the 
selection described in Section 3.2 with the lepton-veto criteria reversed. 



3. The Ur energy spectrum in the b rest frame also depends on whether it is computed 
with HQET or with the spectator model To reflect this difference, the neutrino 
spectrum was distorted with a weight depending linearly on the neutrino energy 
such that the averaged energy changes by ±10%, and the resulting change in the 
measured branching fractions noted. 

4. The Vr energy spectrum depends on the r polarization and, for events containing 
a Ab, on the Ab polarization. These polarizations were varied according to the 
measured uncertainty on the latter, and to the difference between the HQET and 
the spectator model prediction for the former, to determine their effects on the 
branching fractions. 

5. The uncertainty on b — D~X, followed by — >• r^Vj. was propagated to the 
branching fractions. In this process, the missing-energy spectrum is expected to be 
similar to that of the signal. 

6. The small, residual differences between the calibrated missing-energy distributions 
for light-quark events, the b-tagging efficiencies and the lepton-veto effectiveness 
in data and simulation (Section ^) were all attributed to systematic effects and 
translated as such to the branching fraction determination. 
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7. The visible-energy calibration procedure (Section ^4.1| ) is tuned with light-quark 
events, and may not be entirely accurate for bb events in the signal region. In this 
region, the missing energy is found in the simulation to mainly originate from mis- 
reconstructed neutral hadrons. It is observed that bb events yield more hemispheres 
with a large neutral hadronic energy in the data than in the simulation. This excess 
is expected to increase the residual background in the signal region by ~ 20%, 
treated as an additional systematic uncertainty. 

Two additional systematic effects, specific to the b — > r^z/,- D*=^X final state were also 
identified. First, the combinatorial background was estimated in that case from the 
sidebands around the D° peak, and used to derive the corresponding uncertainty on 
BR(b —>■ r^P^D*^X). Second, the b-tagging criterion specific to this channel (Section |3.4| ) 
yields a different efficiency in data and simulation, which affects BR(b r^P^D*^X). 



Table 4: Systematic uncertainties (in %) for BR(b r z^^-X) and BR(b r i'tD*^X). 



Source 



BR(b 



r 



BR(b ^ r-z/^D*±X) 



28 



(xb)=0.702 ± 0.008 
(xc)=0.487 ± 0.008 
BR(b £z/£X)=10.56 ± 0.21 % 
BR(b ^ c ^ iuiX)=7.98 ± 0.22% 



TO. 12 
^0.01 
^0.05 
TO.01 



^0.15 
^0.01 
^0.10 
^0.03 



BR(b ^ D**,D* 



-IT 



=25 ± 7% [|5 



±0.03 



±0.06 



b-fragmentation modelling 



^0.11 



^0.12 



b — ^ r z/^X decay modelling 



±0.06 



±0.06 



BR(D,^ 



=5.79 ± 1.94% 



^0.08 



^0.07 



(P^) = -0.735 ±0.03 
(P(Ab)) = -0.311°:^^ ±0.08^ 



±0.02 
^0.06 



±0.01 
^0.01 



b-tagging efficiency 
/i-identification efficiency 
e-identification efficiency 
Visible-energy calibration 



±0.06 
±0.06 
±0.08 
^0.05 



±0.02 
±0.07 
±0.08 
^0.04 



i?inis« in bb events 



±0.06 



±0.05 



Combinatorial background 



±0.03 



Total systematic uncertainty 



±0.25 



±0.28 



6 Upper limits on B — > t and b — > sviy 

The B^ t^Ut- and b — >■ suu branching ratios are predicted to be too small in the 
standard model to be measured with only four million hadronic Z decays. However, 
upper limits can be set on these branching fractions to constrain possible extensions of 
the standard model, such as type-II two-Higgs-doublet models. Since BR(B^ t^z^t) 
and BR(b suu) can be affected differently depending on the new physics considered, 
the upper limits are conservatively estimated here for each process separately. 
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The missing-energy distribution of the events selected as described in Section ^]3] is 
displayed in Fig. |, and is compared to that of the background. Also indicated in Fig. || 
is the expected enhancement at large missing energy, should either the B~ —>■ T~Ur or the 
b —>■ suu branching ratio be 1%. 




miss 



Figure 6: Missing-energy distributions in the t^Vt and b — > svv final state selection, for the data 

(dots) and for the simulated background (full histogram). Also indicated are the contributions of the 
B~ — > T~ Ur (dashed histogram) and b — s- svD (dotted histogram) processes with a branching ratio of 1%. 
The arrow shows the region in which the two limits are calculated. 



No background subtraction was performed, making the validity of the results presented 
here unaffected by possible systematic uncertainties related to the knowledge of the resi- 
dual background, largely suppressed by a lower cut on the missing energy at 35 GeV. The 
latter cut was optimized so as to maximize the expected 90% CL. upper limit, evaluated 
with simple event counting and in the absence of new physics [^, on the t^v^ 
and b — > sz/i/ branching ratios. The numbers of events, observed in the data and expected 
from background and signal, are displayed in Table |^ in three missing-energy intervals. 
Two events with a missing energy in excess of 35 GeV were observed, with 2.5 ±1.6 events 
expected from all background processes. In absence of any systematic uncertainty, the 
90% CL. upper limits on the t^v^ and b — > si^z/ branching fractions are found to 

be 8.1 X 10~^ and 6.2 x 10^^, respectively. 

However, these limits are affected by the uncertainty on the expected fraction of 
B^ r~z/T- and b — *• ^vv events with such a large missing energy or, almost equivalently, 
with such a large value of Xb. (Most of these events are characterized by a value of Xb 



in excess of 0.9.) This fraction was determined to be 0.1461q;o2i, which translates 
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Table 5: Numbers of events observed in the data in three different £^miss intervals. Also indicated are 
the events expected from the various background and signal processes. For the latter, a branching ratio 
of 1% was assumed. 



-Emiss interval 


[30, 35] GeV 


[35, 40] GeV 


> 40 GeV 


Data 


Selected hemispheres 


31 


1 


1 


Simulation 


Expected hemispheres 


37.0±2.7 


2.5 ± 1.6 


< 1 


Background (detail) 


b t^Ut-X 


7.2 


1.4 




b, c with leptonic decay 


28.9 


1.1 




Other Backgrounds 


0.9 






S 


ignal (BR=1%) 




B^ T^Vj. 


76.9 


53.1 


12.7 


b — >• 


91.0 


78.6 


6.6 



to an uncertainty of 15% on the number of events expected. In addition, the fraction 
of B~ in Z — s> bb events is known [^ to be (38.9 ± 1.3)%. The latter uncertainty of 



3% on the number of events expected affects only the limit on the B~ r^Vr branching 
ratio. These two uncertainties were taken into account following the method of Ref. 
yielding the limits 

BR(B- T-Vr) < 8.3 X 10-^ 
BR(b ^ suu) < 6.4 X 10-^ 

at the 90% confidence level. 



7 Measurement of BR(b r i^t^) with di-leptons 

An alternative method of measuring BR(b —>■ r^P^X) was developed with events where 
both the T and the accompanying D decay to e or /i. Hence, the signature used 
to tag the signal events is a pair of leptons (e, /i) of opposite sign in a jet. The 
background, originating from b ci~U£ followed by c — qi'^i'e, is about 20 times 
larger than the signal. Signal and background are therefore separated on the basis of 
their different kinematic properties. Although this method is statistically less powerful 
than that based on missing energy, it represents an interesting cross-check since it is 
based on a complementary sample of events and sources of systematic uncertainties are 
largely different. Indeed, the main contribution to the systematic error comes from the 
uncertainty on the product BR(b ii^ec, c ii^eq) and on the double charm decay rates, 
B D3D(X) and B ^ DOD(X). 

The presence of three neutrinos in the decay chain of signal events is the main difference 
between signal and background. As a consequence, b — >■ T~UrX decays give larger missing 
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energy, softer lepton spectrum and smaller charged multiplicity to the jet containing the 
lepton candidates. These different kinematic properties of the various categories of events 
are used to separate the signal from the background. A multivariate analysis technique 
with a multilayered neural network (NN) is used to obtain the best discriminating power. 



7.1 Event selection 

Electron and muon identification follows the standard criteria [0, with two refinements 
for muons. Firstly, the muon momentum cut is lowered from 3 GeV/c to 2.2 GeV/c to 
increase the acceptance for the signal (2.2 GeV/c is the minimum momentum for a muon 
to reach the muon chambers) and secondly, any track "shadowed" by another track is 
rejected in the di-lepton selection. Two tracks are said to be shadowing if they share 
in common either hits in the last ten layers of the hadron calorimeter or at least one 
three-dimensional muon chamber hit. This latter cut allows a good description of the 
background to be obtained in the sample of events with two muons in the same jet, as 



shown in Ref. |34|. To obtain the di-lepton sample used for the training of the neural 



network, two additional cuts on the invariant mass Mu of the lepton pair are applied. 

• Electron pair invariant masses are required to be greater than 0.3 GeV/c^ in order 
to reject 7 conversions or Dalitz decays of rj, 7r°. 

• Electron pairs and muon pairs are required to have an invariant mass smaller than 
2.5 GeV/c^ to avoid pairs coming from J/ip decays. 

After these selection cuts, the simulated di-lepton sample consists of 451 signal events, 
7834 h i background events and 5249 di-leptons coming from hadrons mis-identified 
as leptons and light-hadron decays. 



7.2 Analysis method and background estimation 

Jets are reconstructed by the JADE clustering algorithm |p5|| with a i/cut 

of 0.0044 [|17|. 



The following variables are selected as input to the NN: the missing energy in the 
hemisphere of the jet containing the lepton pair, the momenta of the leptons boosted 
to the reconstructed b-hadron rest frame [Q, the invariant mass of the lepton pair, the 



transverse momenta of the leptons with respect to the jet direction, the total energy and 
the charged energy of the jet containing the lepton pair, the fraction of the jet charged 
energy carried by the leptons and the number of charged particles in the jet. 

The main background consists of di-leptons coming from b — ^ c£, c ^ fV^/ decays. To 
determine the product BR(b —>■ i'^'^c) xBR(c —>■ i'lye'), the semi-leptonic branching fraction 
BR(b —>■ ii^eX.) is fixed to the LEP average value |2^, while the h —>■ c —>■ i' fraction is 
determined from the same sample of di-lepton events as follows: a neural network with 
the input variables described above is trained to separate h ci ^ ii'X decays from 
all other processes; and the b — » c — £ fraction is extracted by fitting the simulated 
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neural network output distribution to that of the data. No attempt at an evaluation of 
the systematic uncertainties on this fraction was made. 

Another important background consists of leptons from light-hadron decays or of 
hadrons mis-identified as leptons. A control sample of same sign di-leptons is used to 
test the accuracy of the simulation. This sample has a composition similar to that of the 
background of the analysis, since one lepton candidate is always fake or coming from the 
decay of a light-flavoured hadron. 

The measured values [^] for double charmed b decays (B DsD(X), B D°D(X)) 
are included in the background estimate as well. 



7.3 Branching ratio estimation with a neural network 

The agreement between data and Monte Carlo distributions is good for all the variables 
used in the analysis (Fig. |^). The sample of simulated events is rescaled so that the number 
of selected hadronic events is equal to that in the data. The value of BR(b r~z/^X) is 
obtained by fitting the simulated output neuron distribution to the observed distribution, 
shown in Fig. ^ The discrepancy in the first bin has been traced, using the same 
sign di-lepton sample, to a subset of mis-identified muons not well reproduced by the 
simulation |Q. These background events are well separated from the signal events by the 



NN. The inclusion of the first bin in the fit has a negligible effect on the fitted branching 
fraction. 

The dominant contributions to the systematic errors are the uncertainties in the values 
of the branching ratios involved in the analysis, the uncertainties in the modelling of semi- 
leptonic decays, the effect of the -Emiss calibration previously discussed and those inherent 
to the lepton identification. Some sources of systematic uncertainty (such as those related 
to the lepton identification) which affect both b — > r~i/^X and (b c — * ^) are properly 
taken into account. All the errors were obtained by varying the amount of the component 
of the event sample under study according to its uncertainty, changing the distribution 
of their input variables to the NN when required, and re-fitting the NN output neuron. 
The changes in the fitted fraction of b — >■ T~Ur^ events and in BR(b c ^ i) are 
evaluated with their respective NN, and the corresponding change in the BR(b r^i^rX) 
is calculated. 

The result is 



BR(b ^ r-p^X) = [3.36±0.67 (stat.) ±0.60 (syst.)] %. 



Various consistency checks are performed. The lepton identification performance are 
tested by obtaining results separately for the three subsamples of di-electrons (e''^,e'^), 
di-muons (/i^, /i^) and (e^, /i^) pairs. The measured branching fractions BR(b r'/z^-X) 
are found to be (3.8±1.3)%, (3.6±1.4)% and (2.9±1.0)%, respectively. The simulation of 
the missing energy is tested in two ways. First, a NN fed with charged track information 
only gives BR(b —>■ t^z^t-X) =(3.15 ± 0.89)%. Second, the branching fraction obtained 
when fitting the -Emiss distribution instead of the NN output is (3.05 ±0.90)%. 
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Figure 7: Data and normalized Monte Carlo distributions for (a) the missing energy of the Icptonic 
hemisphere and (b) the momentum of one lepton boosted to the reconstructed b-hadron rest frame. The 
normalization of the simulated distribution for signal events corresponds to the BR(b — > T~r'r^) value 
obtained from the fit. 




0.25 0.5 0.75 1 



Output Neuron 

Figure 8: Output neuron distributions for data and simulation (histograms). The normalization of 
the simulation BR(b t^Vt^) distribution corresponds to the branching fraction value obtained from 
the fit. The samples labelled b T^Vr^ and b /'//fX in the figure contain only the events where 
both leptons are correctly identified. Events with at least one lepton mis-identified are classified as 
"Mis-identification background" . 
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8 Interpretation in two-Higgs-doublet models 



As discussed in Section |I], all five measurements presented in this paper can in principle be 
used to constrain type-II two-Higgs-doublet models. Practically, however, the constraints 
obtained in this framework are largely dominated by two of these measurements. 



• The Z penguin diagrams (Fig. |I|f) are related to photon penguin diagrams already 
severely constrained by b ^ S7 searches, and cannot contribute significantly to 
b [§]. 

• The BR(b r^z/T-D*^X) measurement is significantly correlated with, and is 
statistically less powerful than that of BR(b r^z/^-X). 

• The BR(b — > r~i^T-X) measurement with di-leptons is much less accurate than the 
corresponding measurement based on missing energy. 



Consequently, only the BR(b — > r~z/T-X) measurement based on missing energy and 
the search for the B^ r^v^ final state, kept statistically independent in the analysis, 
were interpreted to constrain type-II two-Higgs-doublet models. In these models, both 
branching ratios depend on the parameter r. 



r = 



tan/3 
mH± ' 



where tan/? is the ratio of the vacuum expectations values of the two Higgs doublets, 
and m^i is the mass of the resulting charged Higgs boson. In the free-quark model, the 
b — *■ r~^',-X branching ratio is expected to be modified with respect to the standard model 
prediction according to p| 



BR(b ^ T-v^X) = BR^^(b ^ r-j/^X) x 



1 - 2m^r^$ + 



mzm 



b^4 



(3) 



where $ is a phase-space factor depending on m,-, rric and amounting to about 0.6. 
An enhancement can therefore be observed for values of r in excess of ~ 0.43 (GeV/c^)~^, 
while the destructive interference yields a reduction of the branching ratio below that 
value. Similarly, the B~ — >■ T~Ur fraction is modified with respect to the standard model 
by 



BR(B 



T Ur, 



brs^(b 



X 



1 2 2 
1 — r m-n. 



(4) 



which shows an even stronger dependence on r, and represents an enhancement for all 
values of r larger than ~ 0.27 (GeV/c^)"^. 

To extract a limit on r, the neutrino energy spectrum and the r polarization were 
computed as a function of r and the simulated events, generated with the standard 
model values, were re- weighted accordingly. The BR(b t^Ut-X.) measurement and the 
BR(B^ t^Vt) limit were derived as presented in Sections ^ and ^, for any r value. 
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From the dependence between BR(b r UrX.) and r of Eq. ^, and from the value of 
T u^X), a 90% CL. upper limit on r was extracted P, ^ 



tan/5/mH± < 0.49 (GeV/c^) ^ at 90% CL. 



(5) 



The T polarization actually plays here an important role since, for r = 0.49 (GeV/c^)~^, 
it amounts to —0.28 while the standard model prediction is —0.735 in b — > r"z/^X decays. 
This effect had been neglected in previous ALEPH analyses |Tl|. 



The dependence between BR(B" 
value of the branching ratio, 



r z/^ 



and r of Eq. H and the standard model 



BRS^(B- ^ T-Ur) = 7.4 X 10^^ (/B/160MeV)^ (|Vub|/0.004)^ , 

were interpreted in terms of a 90% CL. upper limit on r: 

tan/5/mH± < 0.40 (GeV/c^)"^ at 90% CL. (6) 

In the limit extraction, the systematic uncertainties on /b and on Vub p5| were taken 
into account following the method of Ref. I^Sj. The combination of the two results in 
Eqs. ^ and ^ cannot improve on the latter, since the corresponding r value does not 
enhance the b r^z/^X branching ratio. 



9 Conclusion 



With approximately four million hadronic Z decays collected by the ALEPH detector at 
LEP, branching ratios involving a b — > r transition have been measured to be 



BR(b ^ T-jJrX) 
BR(b ^ r-z/^D*±X) 



(2.43 ±0.20 ±0.25)%, 
(0.88 ±0.31 ±0.28)%, 



in agreement with the standard model predictions and consistent with similar measu- 



rements performed by DELPHI |^ and L3 A search for the exclusive decay 

B~ t~Ut- has allowed a 90% CL. upper limit to be set on the corresponding branching 
ratio, 

BR(B- ^ r-z/^) < 8.3 x lO^^ at 90% CL. 

Similar limits have been obtained by CLEO L3 ||] and DELPHI [gg. In the 

framework of type-II two-Higgs-doublet models, these results translate to a constraint on 
the model parameter r = tan/5/mH±, 

tan/5/mH± < 0.40 (GeV/c^)^^ at 90% CL. 

Finally, a limit has been set on the process b suu to 

BR(b sz/z/) < 6.4 X 10"^ at 90% CL. 

A related limit on BR(B^ K^z/z/) has been obtained by CLEO p3[] . 
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